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ABSTRACT 

For an extended sample of Seyfert galaxies we compile from the literature the infrared fluxes in the 
four IRAS bands, the ground-based small-beam (~5-10") fluxes in the standard Q, N, M , L (or L') 
bands, and the nuclear (non-stellar) estimated fluxes in the JHK bands. We estimate nuclear fluxes in 
the L band by applying a correction for stellar light. 

From the statistical study of the infrared colors and luminosities, we derive the typical SEDs of Seyfert 
1 and 2 nuclei and the typical differences in luminosities between the two types of objects in the mid- and 
near-infrared spectral ranges. The analyses of colors and luminosities agree with the fact that in general 
Seyfert 2 nuclei become increasingly fainter with respect to Seyfert 1 as we go toward shorter infrared 
wavelengths. This behavior is consistent with growing anisotropy of a dusty torus emission towards 
shorter wavelengths, but the degree of anisotropy is low (the radiation appears to be substantially 
isotropic at A ^25 /im ). For Seyfert 2 galaxies having Compton-thin obscuring structures at hard 
X-ray energies, we find correlations between the absorbing hydrogen columns and some infrared colors 
and luminosities. 

The observational data appear to severely challenge many models of dusty tori, which hardly account 
for the shapes of the SEDs and the degree of anisotropy observed in Seyfert galaxies. In particular, 
at variance with some earlier claims, very thick and compact tori are basically inconsistent with these 
observational constraints. The most successful models, though having problems in accounting for several 
details, can fit the major infrared observational data of both Seyfert 1 and Seyfert 2 nuclei with tori 
which extend up to several hundreds pc and have fairly low optical thickness. 
Subject headings: galaxies: active - galaxies:nuclei - galaxies: Seyfert - infrared: galaxies 



1 Introduction 

Large amounts of observational data have been accu- 
mulated during the last decade supporting unified models 
of active galactic nuclei (AGNs), particularly of Seyfert 
galaxies. The canonical unified model of Seyfert galax- 
ies claims that Seyfert 1 galaxies (Sis) contain the same 
■central engines as Seyfert 2 galaxies (S2s), but have the 
broad-line region and the strong optical/UV/X-ray con- 
tinuum of the central source concealed from our view by 
a disk or torus of dusty molecular clouds (see, e.g., the re- 
view by Antonucci, 1993). This obscuring structure would 
absorb a significant fraction of the optical/UV/X-ray lu- 
minosity and should therefore radiate this energy strongly 
in the infrared (IR) . This would explain why Seyfert galax- 
ies are strong IR emitters. They appear to be relatively 
brighter than non-Seyfert ones at near-infrared (NIR) and 
mid-infrared (MIR) wavelengths on the basis of the global 
spectral energy distributions (SEDs) and relevant central 
colors (e.g., SpinogHo et al., 1995; Glass & Moorwood, 
1985; Giuricin et al., 1994). 

Dust reprocessing of the primary emission is at present 
regarded as the most likely source of the bulk of the strong 
NIR and MIR emission of radio-quiet quasars and Seyfert 
galaxies. Observational evidence favoring models in which 
the IR emission predominantly or entirely comes from 



heated dust is summarized by Barvainis (1992). 

Geometric, energetic, and spectropolarimetric argu- 
ments indicate that the optical/UV continuum escapes 
anisotropically from the nuclei of at least some Seyfert ob- 
jects. The presence of extended ionization cones which are 
seen in images of emission lines around several S2 galaxies 
(see, e.g., the review by Wilson, 1996) suggests that the 
ionizing radiation emitted from the nucleus is coUimated. 
The observed ionizing photon deficit in S2s (Binette, Fos- 
bury & Parker, 1993) indicates that the central ionized 
gas responds to a higher photon fiux than that seen along 
our line of sight. Broad hydrogen emission lines, presum- 
ably transmitted through the torus, have been detected 
in the IR spectra of some S2 galaxies. (Blanco, Ward & 
Wright, 1990; Goodrich, Veilleux & HiU, 1994; Ruiz, Rieke 
& Schmidt, 1994). These kinds of arguments suggest the 
presence of an obscuring dusty torus in S2 objects (see, e.g., 
the review by Wilson (1992) for a discussion of these top- 
ics). The X-ray spectra of S2s exhibit absorptions larger 
than those of Sis, providing strong support for the pres- 
ence of an obscuring structure (e. g., Nandra & Pounds, 
1994; Smith & Done, 1996). 

The first models of thermal dust emission for the ori- 
gin of the IR emission of Seyfert galaxies and radio-quiet 
quasars simply involve optically thin dust configurations 
(e.g., Barvainis, 1987, 1990, Loska, Szczerba & Czerny, 
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1993), which imply too small an optical depth at opti- 
cal/UV wavelengths (r ^ 3), or configurations in ideal- 
ized slab geometry (Laor & Draine, 1993). Phinney (1989; 
see also Sanders et al., 1989) modeled the IR continua of 
AGNs as thermal emission from a dusty warped disk which 
extends for a few kpc in the host galaxy; but its large ra- 
dial extent would imply an unreasonably large mass, if 
the warped structure were not transparent enough at IR 
wavelengths. 

In recent years radiative transfer calculations specific 
to the IR emission of thick, dusty axisymmetric torus-like 
structures have been carried out by Pier & Krolik (1992, 
1993, hereafter PK), Granato & Danese (1994; hereafter 
GD), Efstathiou & Rowan Robinson (1995; hereafter ER). 
All these models concur that face-on systems (i.e.. Sis 
within the unified scheme of Scyfert galaxies) should be, 
on average, stronger NIR and MIR emitters than edge-on 
systems (i.e., S2s) for a wide range of torus geometries and 
optical depths, whilst at long enough IR wavelengths the 
two types of objects should appear to be similar. However, 
the predicted IR emission and degree of IR anisotropy dif- 
fer in many details. 

As is expected in many models of optically thick tori, the 
observations indicate that Scyfert galaxies display a broad 
IR bump which peaks in the MIR. S2s exhibit steeper IR 
continua between NIR and far-infrared (FIR) wavelengths 
than Sis do (Edelson, Malkan & Rieke, 1987), although 
an accurate determination of the SEDs of the two types 
of objects is hampered at NIR wavelengths by appreciable 
starlight contributions, which are typically greater for S2s 
(Alonso-Herrero, Ward & Kotilainen). Moreover, some 
authors have claimed that S2s are characterized by lower 
mean N-band (A ^10 /im ) luminosities than those of Sis 
(Heckman 1995; Maiolino et al., 1995). The last two facts 
can be interpreted as anisotropy resulting from the pres- 
ence of an obscuring torus. 

Furthermore, the observations show that the MIR spec- 
tra of S2s typically display strong silicate absorption bands 
(at A ~10 yum ), whereas those of SI are mostly featureless 
(Roche et al., 1991). The clusiveness of emission features 
in the MIR spectra of Scyfert galaxies represents a serious 
challenge to many models of tori. 

Comparisons between IR observations and theoreti- 
cal predictions have been presented in several above- 
mentioned papers for fairly small samples of Scyfert galax- 
ies (see also Granato, Danese & Franceschini, 1997). But 
the IR SEDs predicted by torus models and, in particu- 
lar, the level of IR anisotropy deserve to be better tested 
through an extensive statistical analysis of the IR prop- 
erties of Seyfert galaxies. Their NIR and especially MIR 
continua, where the peak of the torus emission is expected 
to occur, are key spectral regions in which to probe the 
structure of the obscuring-reprocessing torus. 

In this paper, in order to test models of IR torus emis- 
sion, we analyze the MIR and NIR photometric data of 
a wide sample of Seyfert galaxies. Specifically, we com- 
pile from the literature IRAS fluxes (at A ~12, ^25, ^60, 
~100 iJ,m ), ground-based small-aperture (~5-10") pho- 



tometric measurements in the standard Q (A ^20 /im ), 
N (A ~10 nm. ), M (A ~4.8 /xm ), L (A ~3.5 /xm ) or L' 
(A ~3.8 fim ) bands, and available non-stellar (nuclear) 
flux estimates in the J (A ^1.25 fim ), H (A ~1.6 fim ), K 
(A ~2.2 /im ) bands. 

In §2 we describe the data sample compiled from the 
literature. In §3 we analyze the MIR, NIR colors and lu- 
minosities of our Seyfert galaxies with the main purpose of 
establishing the relevant SEDs of SI and S2 nuclei and the 
differences in their typical relevant luminosities. We also 
explore relations between these properties and some hard 
X-ray properties. We use survival analysis techniques (see, 
e. g., the review by Feigelson, 1992) in order to exploit the 
information contained in censored data (i.e., upper limits 
on fluxes). In practice we use the software package ASURV 
(Rev 1.2, La Valle, Isobe & Feigelson, 1992), in which the 
methods presented in Feigelson & Nelson (1985) and in 
Isobe, Feigelson & Nelson (1986) are implemented. En- 
suing constraints on the models of IR torus emission arc 
discussed in §4. Lastly, §5 summarizes our results and 
contains our conclusions. 

2 The Data Sample 

We have considered a wide sample of Seyfert galaxies 
(~700 objects) having z < 0.1. Our sample essentially 
contains all known Scyfert galaxies having z < 0.05 and 
many more distant objects. The objects have been taken 
from the lists of Dahari & De Robertis (1988), Whittle 
(1992), de Grijp et al. (1992), Nelson & Whittle (1995), 
and the catalogue of Veron-Cetty & Veron (1996). Like 
several wide samples of Seyfert galaxies considered in the 
literature, our sample is likely to be biased towards objects 
which are bright at optical and infrared wavelengths. 

The Seyfert type classification is primarily based on 
the spectroscopic study of the CfA Scyfert sample by 
Osterbrok & Martel (1993) and on the first four above- 
mentioned papers. The Seyfert types are grouped into the 
three types SI (types 1 and 1.2), SI. 5 (types 1.5) and S2 
(types 1.8, 1.9, 2). The distances D of the nearby galax- 
ies which are included in Tully's (1988) NBG catalog are 
taken directly from the values tabulated therein. The dis- 
tances of non-cluster NBG galaxies have been essentially 
estimated on the basis of velocities, an assumed Ho= 75 
km s~^Mpc~^ and a Virgocentric retardation model in 
which the author assumed that the Milky Way is retarded 
by 300 km from the universal Hubble flow by the mass 
of the Virgo cluster. Non-NBG galaxies are given redshift 
distances D with the same value of Hq. 

We have collected the IRAS fluxes (at A ~12, ~25, ~60, 
~100 /im ) of our Seyfert objects preferentially from spe- 
cific reference sources which report co-added survey data 
(e.g.. Rush, Malkan & Spinoglio (1993), Sanders et al. 
(1995), de Grijp et al. (1992), Edelson et al. (1987), 
Sanders et al. (1989), Mazzarella, Bothun & Boroson 
(1991), Knapp, Bies & van Gorkom (1990), Knapp et al. 
(1989), Soifcr et al. (1989), Bothim et al. (1989), Hill, 
Becklin & Wynn-Wihiams (1988), Rice et al (1988), Helou 
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et al. (1988), Salzer & MacAlpinc (1988), Golombek, 
Miley & Neugebauer (1988)), and from the IRAS Faint 
Source Catalog (see Moshir et al., 1992), in order of pref- 
erence. We have found IRAS data (detections or 3 cr upper 
limits) at A ~12, ~25, ~60, ~100 /xm for 535, 538, 540, 
539 objects, respectively. 

We have considered the recent estimates of the non- 
stellar (nuclear) near- infrared (NIR) J, H, K fluxes. We 
have taken relevant NIR data primarily from the observa- 
tional studies by Kotilainen et al. (1992a, b), who observed 
a hard X-ray-selected sample of Seyfert galaxies (nearly 
all SI), Danese et al. (1992) and Zitelh et al. (1993), who 
observed Markarian SI and SI. 5 objects, Kotilainen & Pri- 
eto (1995), who studied NGC 5252, and Alonso-Herrero, 
Ward & Kotilainen (1996), who investigated S2 objects. 
From NIR imaging observations (often complemented by 
NIR aperture photometry) these authors obtained accu- 
rate estimates of the contribution to the total NIR emis- 
sion arising from the nuclear source alone by dcKiomposing 
the siirface brightness profiles into a combination of a nu- 
clear soiircc, a bulge component and a disc component. 
For NGC 5548 and Mk 509 we have adopted the aver- 
ages of the NIR nuclear fluxes derived by Kotilainen et al. 
(1992a) and Danese et al. (1992). The resulting data sam- 
ple contains nuclear J, H, K fluxes for 73 (41 SI, 12 SI. 5, 
21 S2), 73 (41 SI, 12 SI. 5, 21 S2), and 79 (46 SI, 12 S1.5, 
21 S2) objects respectively. Moreover, we have taken some 
L-band nuclear magnitudes estimated by Alonso-Herrero 
et al. (1996) through model parameters obtained in the K 
band, L aperture photometry and standard K-L colors of 
normal galaxies. 

We have gathered other (less accurate) K-band nuclear 
data from McLeod & Rieke (1995) and from Heisler, De 
Robertis & Nadeau (1996) for 24 and 15 objects respec- 
tively. In their own K-band imaging observations of CfA 
Seyfert objects the former authors have approximated the 
nuclear source as a stellar point spread function scaled 
to the central intensity of the light proflle; this approach 
tends to overestimate the nuclear source contribution. 
Heisler et al. (1996) have reported NIR images of a sam- 
ple of IRAS galaxies whose SED peaks in the IRAS A ~60 
/zm band and have derived the nuclear K-band magnitudes 
of their Seyfert objects from the observed Ha luminosities 
simply by assuming Case B recombination and a power- 
law SED with a spectral index n = —1 {F^ oc i^") in the 
optical and NIR bands. 

In addition, a few JHKL nuclear data have been taken 
from Ward et al. (1987) for NGC 4051 and from McAlary 
& Rieke (1988) for Mk 6, Mk 1040, NGC 2992, NGC 
6814, 1 4329A (for which mean values have been adopted). 
Ward et al. (1987) have fitted the light profile ob- 
tained from R-band CCD images with three components 
(point source, disk, and bulge) and have then extrapo- 
lated the derived nuclear contribution into the infrared 
assuming normal stellar colors. McAlary & Rieke (1988) 
employed a Hubble-law growth curve, modified by a 
color-morphological type-aperture relation, to deconvolve 
multi-aperture optical and NIR photometry of Seyfert ob- 



jects into stellar and non-stellar nuclear emissions. 

Moreover, we have gathered together from the litera- 
ture all published ground-based, small-aperture ('~5"-10") 
photometric measures in the Q, N, M, L or L' standard 
broad bands. L' band measures are considered only for 
objects which lack small-beam L photometry. The consul- 
tation of the " Catalog of Infrared Observations" by Gczari 
et al. (1993) has been of valuable aid in compiling in- 
frared photometric data from the literature. In the case of 
multiple entries for an object, we have adopted the mean 
values of the fluxes published in a reference source. Non- 
detections are reported as 3 cr upper limits on fluxes; in 
case of no detection we have chosen the most stringent 
upper limits on fluxes. 

In general we have preferred the smallest beam size mea- 
sures, in order to minimize the inclusion of starlight in 
MIR fluxes. In any case, we have preferentially selected 
the sets of multi- wavelength observations published by the 
same authors (e.g., the sets of measures by Rieke (1978), 
McAlary et al. (1983), Lawrence et al. (1985), Rudy & 
Rodriguez Espinosa (1985), Edelson et al.(1987). Ward et 
al. (1987), Neugebauer et al. (1987), HiU et al. (1988), 
in order to minimize the effects of data inhomogeneity 
and infrared variability (if any) in the definition of the 
spectral energy distributions (SEDs) of our objects. Long 
time-scale variations of JHKL band fluxes have been fre- 
quently observed in Seyfert galax;ies (e.g., Kotilainen et 
al., 1992b; Glass, 1992), whereas no sure evidence of vari- 
ability at longer infrared wavelengths has been noted in 
Seyfert galaxies and radio quiet quasars (e.g., Edelson & 
Malkan, 1987; Neugebauer et al., 1989). Thus, we have 
paid particular attention to selecting L (or L') fluxes and 
JHK nuclear fluxes measured at the same epoch (if avail- 
able). For this reason we have selected the JHKL nuclear 
fluxes of NGC 2992 and I 4329 A as given by McAlary & 
Rieke (1988) (instead of Kotilainen et al.'s (1982a) values), 
because they are based on small-beam (4.6") JHKL pho- 
tometric observations which are simultaneous (McAlary et 
al., 1983). 

Most of the QNML(L') photometric data have been 
taken from the papers mentioned above and from Stein 
& Weedmann (1976), McAlary, McLaren & Crabtree 
(1979), Rudy, Levan & Rodriguez-Espinosa (1982), Glass 
& Mooorwood (1985), Boisson & Durret (1986), Carico et 
al. (1988), Sanders et al. (1989), Kotilainen et al. (1992), 
Zitein et al. (1993), Wynn-Wilhams & Becklin (1993), 
Vader et al. (1993), Elvis et al. (1994), Maiolino et al. 
(1995), and Alonso-Herrero et al. (1996). We have found 
Q, N, M, L (or L') data (detections or upper limits) for 
45 (20 SI, 5 SI. 5, 20 S2), 173 (65 SI, 17 S1.5, 91 S2), 36 
(14 SI, 7 S1.5, 15 S2), and 202 (79 SI, 21 S1.5, 102 S2) 
objects, respectively. 

If magnitudes are given in the literature, they have been 
generally converted to fluxes by means of Wilson et al.'s 
(1972) zero-magnitude calibration. In view of the fairly 
large photometric errors (i.e., ~10%, ~10%, ~20%, ~20%, 
~25% for the L, L', M, N, Q fluxes), no corrections have 
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been applied for interstellar extinction and redshift. 

Finally, we have considered the measures of 
[OIII]A5007A narrow emission line flux available in the 
literature for 364 objects (105 SI, 35 SI. 5, 224 S2), as 
compiled mostly by Nelson & Whittle (1995), Whittle 
(1992), Dahari & De Robertis (1988), do Grijp et al. 
(1992), Mulchacy, Wilson & Tsvctanov (1996a). 

The starlight contribution dramatically decreases as we 
go from the JHK bands to longer wavelengths (see, e.g., 
the JHKL nuclear fluxes estimated by Ward et al. (1987) 
and McAlary & Rieke (1988)). Also the image sizes of 
luminous infrared galaxies are smaller at A ~3.4 /^m that 
at shorter IR wavelengths (e.g., Zhou, Wynn- Williams & 
Becklin, 1993). The small-beam fluxes in the QNM bands 
of our Seyfert objects can be taken essentially as nuclear 
fluxes. 

On the other hand, the stellar contribution is still ap- 
preciable within small-beam L fluxes; within 3 arcsec it 
ranges between fractions of ^10% to ^^50% (with a mean 
of ~30%) for S2s, which, however, show a greater starlight 
contribution in the JHKL bands than Sis (Alonso-Hcrrero 
et al., 1996). Therefore, wc have corrected the published 
L (or L') fluxes by subtracting an estimated value for 
starlight. This was done by assuming that the J-band flux, 
which has been taken in the same aperture and preferen- 
tially by the same authors as the L (or L') measurements, 
consists entirely of starlight which has J— L=-|-1.3; this is 
the normal color for an old stellar population coming from 
observations of elliptical galaxies and bulge-dominated re- 
gions of spiral galaxies (e.g., Willner et al., 1984). Any 
superposing thermal emission from central hot dust would 
enhance the L-band fltix much more than the J-band flux 
so that the resulting overcorrection in L would be negligi- 
ble. 

Since a few galaxies lack J-band photometry, the re- 
sulting sample of galaxies with L-band nuclear fluxes or 
relevant upper limits comprises 200 objects (78 SI, 21 S2, 
101 S2). We have checked that our rough method would 
give L nuclear magnitudes similar to those less crudely de- 
rived by Alonso-Herrero et al. (1996) for 11 galaxies (in 
this case we directly took their values). In a few cases our 
L-band correction is slightly greater than the measured L- 
band fluxes. For these galaxies, which have probably very 
little or no L-band flux associated with an active nucleus, 
we have adopted an upper limit of F^^ < 1 mJy. 

The large-beam IRAS MIR fluxes at ~12 and ~25 
yum certainly include emission coming from the host galax- 
ies. However, their good correlations with the K nuclear 
emission and the N emission suggest that the bulk of 
IRAS MIR emission of Seyfert galaxies comes from the 
nuclear region, at variance with the IRAS FIR emission 
(at A ~60 and ~100 fim ) which mostly originates from 
the host galaxies (see, e.g., Danese et al., 1992; Giuricin, 
Mardirossian & Mezzetti, 1995). This is also directly re- 
vealed by the high-resolution IR images of the central re- 
gions of nearby Seyfert galaxies, since these observations 
clarify the relative role of AGN and star formation in pro- 
ducing the IR emission (see e.g. the review by Telesco 



(1988) and the MIR maps of NGC 1068 (Telesco & Decher, 
1988), NGC 1808, NGG 1365, and NGC 5506 (Telesco, 
Dressel & Wolstencroft, 1993), NGC 7469 (Miles, Houck 
& Hayward, 1994)). Therefore, in this paper we mainly 
deal with the IRAS MIR fluxes. We use the IRAS FIR 
fluxes only for some comparisons. 

Throughout this paper the fluxes measured in the IRAS 
A ~25 and ^^12 i^m , Q, N, M bands and the non-stellar 
(nuclear) estimated fluxes in the L, K, H, J bands are re- 
spectively denoted by F25, F12, Fq, Fn, Fm, Flu, Fku, 
Fhu, Fjn- The corresponding monochromatic luminosi- 
ties Lj/ (expressed in W/Hz) and the luminosity L^qiii] 
(in erg/s) of the [OIII]A5007 

A emission line follow from the adopted distances (with 
Ho=75 km s~^ Mpc^^). We have always used L'-band 
measures in place of L-band measures, if the latter are not 
available. We have verifled that SI objects alone give re- 
sults (color and luminosity distributions) similar to those 
of Sl-t-Sl.5 objects. Therefore, in the following sections 
the latter objects will be simply denoted as Sis. 

3 Analysis and Results 
3.1 The IR colors 

In order to establish the behavior of the IR 
continuum at MIR and NIR wavelengths, we 
have evaluated the following colors (decimal log- 
arithm of flux ratios) log(F25/fi2), log(FQ/FAr), 

l0g{FN/FM), \0g{FM/FLn), log(f jv/Fl„), log(Fjv/F^„) , 

log{FLn/FKn), iog{FKn/ Fhu) , ^og{FHn/Fjn) Separately 
for S1+S1.5 and S2 galaxies. We have evaluated the nine 
above-mentioned colors (or relevant lower or upper limits) 
for all objects detected in the two relevant bands and also 
for the objects having upper limits on F12, Fq, Fm, Fm, 
Fn, Fn, for the flrst six colors, respectively. The very few 
objects having upper limits on the other term of a given 
ratio or those which have upper limits on both terms are 
excluded. There are only detections for the case of the 
last three colors. 

We have employed the Kaplan-Meier product limit es- 
timator (which is part of the ASURV package) in order 
to calculate the median, mean and distribution functions 
of the colors of the two groups of Seyfert galaxies. Re- 
distributing upper limits (if any) uniformly along all the 
intervals of lower detected values, the Kaplan-Meier esti- 
mator provides an efficient, non-parametric reconstruction 
of information lost by censoring in the case of randomly 
censored data sets. We compare the distribution func- 
tions of colors for the two classes of Seyfert galaxies by 
testing the "null" hypothesis that the two independent 
random samples are randomly drawn from the same par- 
ent population. In the case of uncensored data, this is 
accomplished by using the Kolmogorov-Smirnov test (see 
e.g. Hoel, 1971). In the general case of censoring we use 
two versions of Gehan's test (one for permutation variance 
and the other with hypergeometric variance), the logrank, 
the Peto-Peto, and Peto-Prentice tests. 
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Fig. 1. — The histograms show the color distributions for SI (solid lines) and S2 nuclei (dashed lines). 
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These two sample tests give; the probability p that the 
two data samples would be drawn from one parent popula- 
tion. The statistical significance of the difference between 
the color distributions of the two classes of galaxies is at 
the 100 (1 — p)% level. The two-sample tests differ in 
how the censored points are weighted and consequently 
have different sensitivities and efficiencies with different 
distributions and censoring patterns. 



3.1.1 The color distributions 

Table 1 gives the numerical outcomes, namely the total 
number of objects N, the number N^i of relevant upper 
limits, the median and mean values (with associated 1 a 
errors), and the corresponding median spectral index n 
(Fy cx z/~") for the distributions of the above-mentioned 
nine colors. Fig. 1 shows the Kaplan-Meier distributions 
of the colors for SI and S2 objects. 



In general S2s have redder log(F25/i^i2), logiFq/FN), 

\0g{FN/FM), \0g{FM/FLn), log(i^jv/i^L„), log(i^jv/i^/f„), 

\og{FLn/ Fku) colors (i.e., steeper SEDs) than Sis. Table 
2 contains the mean and median values of the probability 
p that the colors of SI and S2 objects would be drawn from 
the same parent population. One can see that the statis- 
tical significance of the color differences between the two 
groups of galaxies is high in several cases and marginal in 
the cases of log(^Q/i^Ar) and log(f at/Fm). On the other 
hand, the two groups of Seyfert objects do not differ ap- 
preciably in the log (Frth /-?>/„) and \og{FHn/Fjn) distri- 
butions. 

We have attempted to investigate the possible effects of 
circumnuclear star formation on IR colors. Although it is 
difficult to come up with consistent quantitative indica- 
tors of this phenomenon, it does not seem to occur very 
frequently in nearby Seyfert objects (e.g., Pogge, 1989). 
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TABLE 1 
The IR Colors 



Color 


Sample 


N 
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0.37± 0.04 


1.4 




82 


12 


1 


^ -™-i-n ofi 

— 0.15 


0.50± 0.07 


1.8 


log(FM/-Fi„) 


81 


23 


2 


„ ^ n + O 02 


0.18± 0.03 


1.3 




82 


16 





„ nrr-l-O 13 


0.39± 0.07 


2.0 


log{FN/FKn) 


81 


38 


4 


„ „„-i-0 03 


0.85± 0.04 


1.3 




82 


16 


1 


1 /-n + 0.18 

l-69Io.53 


1.56± 0.16 


2.6 


log(Fi„/FKn) 


81 


30 





36+° °'' 


0.32± 0.05 


1.8 




82 


11 





^■'°-0.10 


0.74± 0.11 


3.7 


log(FK„/FH„) 


81 


57 





r, 97+0.02 
"■^'-0.01 


0.28± 0.03 


2.0 




82 


22 





19+0-15 


0.25± 0.05 


1.4 




SI 


57 







0.22+ 0.02 


1.9 




S2 


22 





''■^■'-0.02 


0.22± 0.08 


2.1 



In any case, the galaxies which are known to have clear 
signs of intense circumnuclear star formation (e.g., NGC 
1068, NGC 1365, NCC 1808, NCC 3982, NGC 4388, NGC 
4945, NGC 5427, NGC 5728, NGC 5953, NGC 6221, NGC 
7469, NGC 7582, NGC 7592, I 4553, Circinus galaxy), as 
is evidenced mostly by spectroscopic investigations of nu- 
clear spectra or emission line imaging surveys (e.g. Veron- 
Cetty & Veron, 1986; Pogge, 1989, Evans et al., 1996), do 
not appear to have colors significantly different from the 
norm. 

TABLE 2 

The p-values for the IR Color Distributions 



Color 


Mean p 


Median p 


(1) 


(2) 


(3) 


log(F25/Fi2) 


< 0.0001 


< 0.0001 


log{FQ/FN) 


0.138 


0.087 


\og{FN/FLn) 


< 0.0001 


< 0.0001 


logiF^/ Fm) 


0.059 


0.065 


log{FM/FL„) 


0.007 


0.009 


log(Fjv/FKn) 


0.002 


0.0001 


log(FLn/FK„) 


0.0009 


0.0009 


log(Fx„/Fjy„) 


0.21 


0.21 


log(F/f„/Fj„) 


0.63 


0.63 



The color distributions are in general Gaussian and uni- 
modal, with some exceptions mentioned below. We have 
adopted the probability Pw associated with the W-test 
(Shapiro & Wilk, 1965) in order to estimate the Gaussian- 
ity of the color distributions for uncensored data. 

We have computed the dispersions of the color distri- 
butions for the uncensored data only. We have compared 
these observational dispersions aobs with the values acai 



which arc expected from reasonable observational uncer- 
tainties on the fluxes, using the first order approximation 
a{logx) ~ log{e)a{x)/x. We have adopted the following 
percentage errors: 10% on F25, 15% on F12, 25% on Fq, 
20% on Fa,, 20% on Fm, 15% on 30% on Fku, Fhu, 

Fjn- 

As shown by Table 3, the expected values {<Jcai) are al- 
most always smaller than the observed ones {(Jobs)- This 
indicates the presence of appreciable intrinsic dispersion 
in the IR colors. The observed dispersions for S2s are typ- 
ically greater than those for Sis even at the longest IR 
wavelengths, where uncertainties in the fluxes are similar 
for the two types of objects. 

Some comments on the individual color distributions fol- 
low. 

3.1.2 The mid-IR colors 

Let us discuss the MIR colors log(F25/Fi2) and 
logiFq/ Fn). The distribution of log(F25/i^i2) for Sis 
presents a tail at low values which makes its shape a little 
different from a Gaussian one (the W-test applied to the 
uncensored data gives a probability Pw=0-026 that the 
distribution does not deviate from Gaussianity) . For Sis, 
our log(F25/-F'i2)-distribution is shifted to somewhat lower 
values than previous distributions based on much smaller 
data samples, whilst it essentially confirms previous typ- 
ical results for S2s (Giuricin et al., 1995; Maiolino et al., 
1995). 

The log(FQ/FAf)-distribution, as yet unexplored in the 
literature, has mean and median values similar to those of 
the log(F25/Fi2)- distribution for S2s. On the other hand. 
Sis lead to a log(F(3/F]v)-distribution which is shifted to 
somewhat greater values than the log(F25/Fi2) one. 

TABLE 3 

The Dispersions of the IR Color Distributions 



Color 


<^obs 


<^obs 






81 


82 




(1) 


(2) 


(3) 


(4) 


log{F25/Fi2) 


0.19 


0.21 


0.08 


log{FQ/FN) 


0.12 


0.25 


0.14 


iog{FN/FLn) 


0.18 


0.40 


0.11 


\og{FN/FM) 


0.17 


0.25 


0.12 


Iog(FM/Fi,„) 


0.16 


0.27 


0.11 


log{FN/FKn) 


0.23 


0.75 


0.16 


log(Fi„/FK„) 


0.27 


0.37 


0.15 


^Og{FKn/FHn) 


0.19 


0.25 


0.18 


\og{FHn/Fjn) 


0.13 


0.39 


0.18 



The evaluation of the log(F25/Fi2) distributions for the 
25 Sis and 17 S2 objects with available log(FQ/Fjv) con- 
firms the presence of some systematic differences for Sis, 
since their distribution leads to a median value of 0.31 
and a mean value of 0.32. Consistently, a steepening of 
the MIR SED of SI objects towards shorter wavelengths 
is suggested by Roche et al.'s (1991) MIR spcctrophoto- 
metric observations. These authors noted that the slopes 
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of the A ^^8-13 /im spectra they obtained on ground-based 
telescopes with small beams (~ 5" ) were steeper than the 
IRAS 12-25 iim slopes for galaxies with featureless MIR 
spectra (mostly Sis). 

3.1.3 The mid-to-near IR colors. 

The seeming bimodality of the log(_FAr/FM)-histogram 
is not statistically significant. The \og{FM / Flv) dis- 
tribution for S2s presents a tail at high values, which 
makes it significantly non-Gaussian. (Pn/=0.002). The 
\og{Ffq / Fm) and \og{FM / Flu) distributions give results 
which are substantially consistent with those relative to 
the log(F/v/FL„) distributions, since they imply median 
(mean) values of log(F/v/i^i„) equal to 0.63 (0.55) and 
0.84 (0.89) for Sis and S2 objects, respectively. The 
log(FAr/_F'L„)-distribution for S2s is significantly more 
peaked than a Gaussian distribution {Pw < 0.0001). In 
the following we shall adopt the log(FAr/FL„)- distribu- 
tion in order to define the SED between the relevant wave- 
lengths, because the statistics arc better and there is no 
significant indication of a change (at A ~4.8 /zm ) in the 
SED slopes of the two types of objects. Wc point out that 
neglecting a correction for starlight in the L-band, one 
obtains an appreciable overestimate of this color also for 
small-beam observations of nearby objects. For instance, 
the log{Fi^/FL) distributions of 67 SI and 66 S2 objects 
lead to median (mean) values of 0.44 (0.42) and 0.74 (0.72), 
respectively. Compared to the results of the corresponding 
log(i^Ar/-FL„)-distributions, this implies a median L-band 
starlight contribution of 28% and 55% for the two groups 
of objects, respectively. Correction for starlight is gener- 
ally neglected in the relevant literature, except in Pier & 
Krolik (1993). Encouragingly, relying on the small data 
sample (15 Sis and 6 S2s) considered by these authors, we 
obtain similar medians of log(F/v/FL„) (i. e., ~ 0.6 and 
0.9, respectively). 

The two fluxes used to compute the quantities 

l0g(F25/i^l2), \0g{FQ/FN), \0giFN/FM), \0g{FM/FLn), 

log(i^Ar/i^i,„) come from the same reference sources in 
many cases. This limits the effect of data inhomogene- 
ity on these five colors. This problem may be more se- 
rious for the sample of \og{FM / Fku), in which the two 
fluxes generally come from different reference sources and 
refer to measurements carried out in different epochs. We 
have obtained a log{F]\[ / F^n) distribution characterized 
by median (mean) values of 0.87 (0.87) for 51 Sis (with 6 
upper limits on Fn) and 1.06 (1.30) for 38 S2s (with 4 up- 
per limits on Fn). But wc have preferred to rely finally on 
a subsamplc of bona fide \og{FM / Fku) colors, by consid- 
ering only objects with Fku data taken from the accurate 
decomposition analyses by Kotilaincn et al. (1992a, b), 
Danese et al. (1992), Zitelh et al. (1993), Kotilainen & 
Prieto (1995), and Alonso-Herrero et al. (1996). Table 1 
and Fig. 1 present our results for this selected subsample, 
which, compared to the above-mentioned wider sample, 
substantially confirms the median and mean values of the 
\og{Fjq / Fxn) distribution for Sis, but leads to redder col- 
ors for S2s. The addition of a few objects having Fku 



estimates taken from Ward et al. (1987) and McAlary & 
Rieke (1988) yields a \og{FN / Fku) distribution which is 
similar to that of the selected subsample. On the other 
hand, the inclusion of Fku fluxes taken from McLeod & 
Rieke (1995) and Heisler et al. (1996) tends to shift the 
log{FN / Fku) distribution of S2s to appreciably smaller 
and larger values, respectively, and gives rise to an even 
broader and more irregular (i.e., non-Gaussian) color dis- 
tribution than the broad one of the selected subsample. 
The large proportion of Fku data taken from McLeod & 
Rieke (1995) in the wide sample accounts for the above- 
mentioned difference in the median colors of S2s. 

3.1.4 The near-IR colors 

The sample oi \og{FLn/ Fku) colors leads to a distribu- 
tion characterized by median (mean) values of 0.34 (0.34) 
for 66 Sis (with no censoring) and 0.65 (0.75) for 44 S2s 
(with one upper limit on -Fl„). But this sample suffers 
from the same problems of inhomogeneity which affect 
the log{FN / Fku) colors. Therefore, also in this case we 
have preferred to rely on a selected subsample of bona 
fide log{FLn/ FKn) colors, by considering only objects with 
both good Fku data (defined as above) and with F^n 
fluxes derived from L band data measured by the same 
authors in similar epochs; the latter requirement is in- 
tended to avoid the possible effects of nuclear variability. 
Table 1 and Fig. 1 refer to this selected subsample, which, 
however, substantially confirms the typical values given by 
the wider sample for both Sis and S2s. These two distri- 
butions were found to be not significantly different from 
those of the selected subsamples through the application of 
the usual two-sample tests, although the selected subset of 
S2s seems to lead to a narrower spread of values and that 
of Sis shows a tail at low values (Ph'=0.007). Incidentally, 
a somewhat larger subsample which contains only objects 
with good Fku data (as above), but with no restrictions 
on the reference sources for the F^n data, again essentially 
confirms the typical values listed in Table 1. 

The quantities log(Jjv/Ji,„)- and log{FLn/FKn) yield 
results which are in substantial agreement with those given 
by the log(i^Ar/PVra)-distribution, since for this quantity 
they imply median values of 0.91 (for Sis) and 1.68 (for 
S2s). The SEDs of both types of objects considerably 
steepen beyond the L-band (i.e., at shorter wavelengths). 
In the following we shall use the colors log(Jjv/i^_L„)- 
and log{FLn/FKn) to establish the SEDs between relevant 
wavelengths. 

As regards the quantities ^og{FKnl Fhu) and 
\og{FHn/ Fjn), S2s show distributions which are broader 
than those of Sis, but they are centered at similar values 
and have no significant bimodality. We have verified that 
the distributions of theses two quantities do not change ap- 
preciably, if we restrict these two data samples to objects 
with Fxn data taken from the accurate decomposition 
analyses by Kotilainen et al. (1992a, b), Danese et al. 
(1992), Zitelh et al. (1993), Kotilainen & Prieto (1995), 
and Alonso-Herrero et al. (1996). 
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TABLE 4 
The IR Luminosity Distributions 



13.5 



log V 



14 
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Fig. 2. — This figure shows the median SEDs of SI (solid 

lines) and S2 nuclei (dashed lines). The two SEDs are ar- 
bitrarily normalized so that they match at long wavelengths 
(A = 20 yum , Q-band). 

The absence of any statistical difference between Sis 
and S2s in these two colors is clearly unexpected. It is 
likely to reflect very large uncertainties associated with 
the delicate estimates of nuclear fluxes in the H and J 
bands for S2s, where the starlight contribution is even 
greater than in the K band (see, e. g., Alonso-Herrero 
et al., 1996; Kotilainen et al., 1992a). Therefore, in the 
following we will adopt take into consideration only the 
results relative to Sis, disregarding those relative to S2s. 

Fig. 2 illustrates the adopted median SEDs of SI and S2 
nuclei as they result from the colors discussed above. The 
two SEDs are arbitrarily normalized so that they match 
at long wavelengths (A=20 /xm ). Since S2s always have 
steeper SEDs than the former objects, the degree of dif- 
ference between the luminosities L^, of the two groups of 
objects increases towards shorter wavelengths. 

3.1.5 Correlations involving colors 

We have investigated the correlations between colors by 
computing the two non-parametric rank correlation coefS- 
cients, Spearman's rg and Kendall's rk (see, e.g., Kendall 
& Stuart, 1977), in cases of uncensored quantities. In 
testing the significance of the correlations between the 
variables x and y containing censored data, we have re- 
lied on the Cox proportional hazard model, the general- 
ized Kendall rank correlation statistics, and the general- 
ized Spearman rank-order correlation coefficient. The first 
method allows censored data for only one variable, whereas 
the last two methods also allow simultaneous censoring in 
both variables. 

As expected, the quantity log(F25/Fi2) appears to cor- 
relate significantly with \og{FQ/ F^) for both Sis and S2s 
(at the median, one-tailed significant levels of 96.5% and 
99.2%, respectively). 
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SI 
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S2 


91 


14 




22.58± 0.08 
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SI 
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SI 
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22.65± 0.07 




S2 
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8 
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'■-'-O.OB 


21.84± 0.08 


log Lku 


SI 


58 





22 oq+0.19 


22.41± 0.07 




S2 


21 





2j 10+0.21 


21.09± 0.17 



In general we have detected no other significant correla- 
tions, with the exception of a log(F25/.F'i2)- log(-Fiv/Fi„) 
correlation for SI galaxies (N=68, 99.0% significance 
level), which holds also for uncensored data. 

We have verified that no color correlates with the galaxy 
distances (for both Sis and S2s). Subsamples of nearby 
objects give similar color distributions. 

Moreover, we have searched for correlations between col- 
ors and some indicators of the MIR luminosities such as 
Ln, L12, L25, because MIR emission, being an approxi- 
mate constant fraction of the bolometric flux of Seyfert 
nuclei, is believed to provide the best indication of the 
nuclear bolometric luminosity (Spinoglio & Malkan, 1989; 
Spinoglio et al., 1995). In general, we have detected no cor- 
relations, except for the color log(i^25/-F'i2); which shows 
a slight tendency to be somewhat greater in SI and S2 ob- 
jects of greater MIR luminosities. This behavior is dissim- 
ilar from that shown by log(F6o/F25) and log(i^ioo/-F6o)i 
which tend to decrease with increasing MIR luminosity in 
Seyfert galaxies. The latter behavior could be ascribed 
to a decreasing relative contribution of the cool emission 
from galaxy disks in high-luminosity AGNs (Giuricin et 
al., 1995). 

3.2 The IR luminosities 

We are interested in establishing the differences be- 
tween the typical luminosities of Sis and S2s in differ- 
ent bands, because appreciable differences are predicted 
by torus models (especially for very optically thick con- 
figurations and/or at the shortest IR wavelengths). For 
the sake of comparison with samples published in the lit- 
erature, in Table 4 we list the total number of objects N, 
the number of upper limits Nui , the mean and median val- 
ues of Li, together with associated la errors, for the two 
groups of objects and different spectral bands. 

Table 4 reveals that Sis tend to be more luminous than 
S2s at the shortest IR wavelengths, but this may be par- 
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tially due to observational selection effects which favor the 
observation of luminous objects at great distances (Sis are 
typically more distant than S2s). In order to remove this 
distance selection effect in the comparison of groups of ob- 
jects lying at different typical distances, for each spectral 
band we have considered the logarithmic relation 



logL^ = 17.078 + log{F^,min) + 2logD 



(1) 



which yields the minimum value of logL (with L,^ in 

W/Hz) corresponding to the minimum uncensored or cen- 
sored value of (-Fi/,mm) for an object at distance D 
(in Mpc). Obviously, upper limits and the detections 
of the faintest objects cluster near this relation in the 
log Li, — log D diagram. In the case of Lum we have taken 
the subset of objects with Fku data taken from the accu- 
rate decomposition analyses by Kotilainen ct al. (1992a, 
b), Danese et al. (1992), Zitelli et al. (1993), Kotilainen 
& Prieto (1995), and Alonso-Herrcro et al. (1996). 

We have used relation (1) to evaluate the observedmiims 
calculated vahic of logi^, i. e. A(log_L,y), for each object. 
In this manner we have obtained A(log Ly)-distributions 
that are free from distance selection effects. Table 5 lists 
the total number of objects N, the number of upper lim- 
its Nui, the mean and median values of A{\ogL^) (to- 
gether with associated 1 a errors), for the two groups 
of objects and different spectral bands. Fig. 3 dis- 
plays the A(log Lix)-distributions for the most interesting 
cases. There is no significant evidence of bimodality in the 
A(logL^)-distributions. The statistical significance of the 



difference between the distributions of the two groups of 
objects has been evaluated through the application of the 
above-mentioned five two-sample tests. In several cases 
the non-random character of the censoring pattern tends 
to give asymmetrical distributions. However, since two- 
sample tests seek only to compare two samples, they do 
not require that the censoring patterns of the two samples 
be random. 

Table 5 and Fig. 3 reveal that S2 galaxies have, on aver- 
age, greater L25 luminosities than Sis by a median factor 
of~1.5 (at a significance level corresponding top <0.0001). 
The Lq distributions give consistent results, although the 

difference is never statistically significant, because of the 
relatively low number of objects. As we go towards shorter 
wavelengths the difference between the typical luminosities 
of the two classes first disappears (see, e.g., L12 and L^v) 
and then reappears in the opposite sense (see L^n and 
Lku)- As a matter of fact, the Llu luminosities of S2s 
are significantly fainter than those of Sis by a factor 1.8 
(with a median p=0.0004). In other words, as we go from 
^25 to Lin, the ratio of the S2 luminosities over SI lumi- 
nosities typically decreases by a factor of ~2.7. The Lku 
luminosities of S2s are signific;antly fainter than those of 
Sis by a factor ^3.5 (with a nicdian p <0.0001); the total 
sample of Lku data gives ^3.7. Hence, in going from L25 
to Lku, S2s become typically fainter than Sis by a factor 
~5.2 (~5.5 for the total sample). We have verified that 
subsamples of nearby objects always give similar results 
for all spectral bands. 
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Fig. 3. — The histograms show the distributions of A(logL,^) relative to different photometric bands, for SI (solid lines) and S2s 
(dashed lines). A(logL„) is the observed minus calculated value of the logarithm of the monochromatic luminosity. 
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TABLE 5 

The A{logLi,) Luminosity Distributions 



Variable 


Sample 


TV 


N^i 


Median 


Mean 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A log L25 


SI 


203 


45 


r, 70 + 0.04 

^■'=-0.04 


O.Slit 0.04 


A log L12 


S2 
SI 


335 
203 


60 
69 


0-951^3 
0.82l»:°3 


l.Olit 0.03 
0.86± 0.04 


A log Lq 


S2 
SI 


332 
25 


131 



0.82l^:« 


0.87± 0.04 
0.61± 0.07 


A log Ljv 


S2 

SI 


20 

82 


3 

15 


n 77+0-21 

' -0.09 

n /AO + 13 

0-^S-om 


0.85± 0.14 

0.97± 0.06 


A log Lm 


S2 
SI 


91 
21 


14 



90+"-" 

^•"■'-0.03 


0.96± 0.07 
0.60± 0.08 




S2 


15 





Q 70+0.10 
'J- ''^-0.15 


0.82± 0.15 


A log Ll„ 


SI 


99 


2 


1 cc+0.07 


1.66± 0.06 


A log Lkh 


S2 
SI 


101 
58 


8 



-^•^'^-0.10 
9 ,7+0.09 
^•^'-0.10 


1.37± 0.07 
2.10± 0.06 




S2 


21 





1 62+°" 


1.64± 0.16 



Interestingly, the increasing dimming of S2s towards 
shorter wavelengths (with respect to Sis) agrees with 
the differences between the SEDs of Sis and S2s and 
with the expectations of models of obscuring tori emit- 
ting anisotropic radiation at short enough wavelengths. 
But it is a little surprising that S2s appear to be typically 
brighter than Sis at A ~25 /xm , because in the frame- 
work of torus models S2s are expected to have similar or 
smaller luminosities. By the way, in the same manner we 
have found that S2s arc typically brighter than Sis also 
at IRAS A ~60 and ~100 /im bands by similar median 
factors, i.e ^1.5 and ~1.7, respectively. 

We have deemed it useful to check this point by con- 
sidering some subsamples, namely the UV excess-selected 
Markarian objects, the IRAS 12 /zm flux-limited Seyfert 
galaxies listed by Rush et al. (1993), and the IRAS sam- 
ple by de Grijp et al. (1992), who selected objects with 
relatively warm 25 /xm to 60 /um colors. Only the first and 
second subsamples confirm the dissimilarity between the 
A(logZ/y)-distributions of Sis and S2s in the three above- 
mentioned IRAS bands. This dissimilarity turns out to be 
fairly strong in the first subsample (where it amounts to 
a factor of ~1.6 for A(logL25)), weak in the second sub- 
sample (where the corresponding factor is ~1.3), whilst it 
disappears in the third one. Hence, there is no sure evi- 
dence of an intrinsic difi^crcncc between the two classes of 
objects. Our results may simply reflect a selection effect 
in several current samples, in the sense that S2s under- 
going active star formation are more likely to be detected 
at short IR wavelengths and especially at optical wave- 
lengths than objects not undergoing star formation. Tims 
star formation may seem more common in galaxy disks 
hosting S2 nuclei than in those containing Sis, as recently 
stressed by Maiolino et al. (1995), who examined the 
extended MIR emission of Seyfert galaxies by comparing 
the A ~10 ^m ground-based observations with the IRAS 
A ^12 /im fluxes. The aforementioned selection effect may 
also explain why circumnuclear star formation seems to be 



a little more common in S2s than in Sis, as suggested by 
the emission line imaging surveys of nearby Seyfert galax;- 
ies (e.g., Pogge (1989) and Evans et al. (1996)). 

We have provided a further check on this point by com- 
paring the distributions of the IR luminosities normalized 
to a quantity which is powered by the active nucleus and is 
presumed to be isotropic, i.e free of viewing angle effects, 
like optical emission lines produced in the NLR, hard X- 
ray emission and central radio continuum emission (e.g., 
Mulchaey et al., 1994). We have chosen to use the lu- 
minosity of the [OIII]A5007 A emission line, because it 
gives the widest data sample, although this quantity can 
be affected by NLR dust absorption (e.g., Dahari & De 
Robertis, 1988a) and torus obscuration, as was suggested 
in the case of radio galaxies (Jackson & Brown (1990), Hes, 
Barthel & Fosbury (1993), Laing et al. (1994)). Although 
extended line emission may be missed in small-aperture 
measurements, these measures are less dependent than ra- 
dio fluxes on the size of the observed region (the central 
radio emission of Seyfert galaxies has been recently dis- 
cussed by Giuricin, Fadda & Mezzetti, 1996). 

Sis and S2s have log(L[o/77])-distributions character- 
ized by mean (median) values of 41.15 (41.22) for N=140 
objects and 40.90 (41.02) for N==224 objects (with L^oiii] 
in erg/s), together with A log(_L[Q///])-distributions char- 
acterized by corresponding values of 1.55 (1-57) and 1.56 
(1.52). The two Alog(L[o///])-distributions do not dif- 
fer, which is in line with the common view that the 
two types of Seyfert galaxies have the same [OIII] prop- 
erties (e.g., Dahari & De Robertis, 1988a; Keel et al., 
1994; Mulchaey et al., 1994; Mulchaey, Wilson & Tsve- 
tanov, 1996b). Comparing the distributions of the quan- 
tities l0g(Li2/L[O/77]), l0g(iQ/i[O777]), log(i25/i[07//]), 

log(L6o/i[0777]), log(iioo/i[0777]) by means of the usual 
two-sample tests we have again found that S2s have 
brighter normalized luminosities than Sis, in the last 
flve bands, by similar median factors (i.e. ~1.5, ~1.4, 
'^l.B, '^1.7, respectively), whereas they do not differ 
from Sis at A '--^12 /zm . The analysis of the quanti- 
ties l0g(L]v/L[0777]), l0g(-E'Ln/-^[0777]), log(ijfn/.i^[0777] ) 

shows that S2s have fainter normalized luminosities than 
Sis by median factors of ~1.5, ~2.4, ~5.9, respectively. 
Again we have verified that subsamples of nearby objects 
give similar results. 

Incidentally, within extended data samples (though 
smaller than ours) examined in the literature, de Grijp et 
al. (1992) and Keel et al. (1994) found similar FIR lumi- 
nosities for the SI and S2 objects of their IR-selected sam- 
ple, whereas Dahari & De Robertis (1988) and Mulchaey et 
al. (1994) noted enhanced normalized IRAS luminosities 
(at least in one of the above-mentioned three bands) for 
the S2s of their samples which are dominated by Markar- 
ian objects. All these results are consistent with what 
we found for the overall sample and subsamples based on 
different selection criteria. 
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TABLE 6 

Correlations between Njj and IR Properties 



m 1.0 
o 

I 

22 23 24 
log (cm"^) 

Fig. 4. — The plots show the correlations between in- 
frared quantities and the absorbing hydrogen colums Nh (ex- 
pressed in cm~^). The infrared luminosities are integrated 
luminosities evaluated as u L„. 

In any case, our results clearly favor models in which 
obscuring structures are essentially isotropic emitters al- 
ready at A ~25 m. It is encouraging that the degree of 
the increasing faintness of S2s towards shorter wavelengths 
(with respect to Sis) - as derived from a comparison of the 
luminosities - turns out to be in fairly good, quantitative 
agreement with what we expect from the differences be- 
tween the SEDs of the two groups of objects. S2s are 
found to be fainter than Sis by a typical factor of ~1.5 
at A ~12 iim , of ~2.7 at A ~3.5 /xm , of ~5-7 at A ~2.2 
/xm . (Specifically in the last case, values of ~5.2 and ~6.6 
are derived from the differences in the luminosities and in 
the SEDs, respectively.) The analysis of the IR luminosi- 
ties normalizcid to L[o///] analysis which wo regard 
as less reliable than the previous ones - gives similar re- 
sults at long IR wavelengths and gives a somewhat more 
pronounced dimming of S2s at short IR wavelengths. 

Within the framework of torus models these numbers 
establish the degree of anisotropy of the IR emission. 
Notably, we find a lower degree of anisotropy than that 
predicted by many torus models (see next section) and 
claimed by Heckman (1995). For instance, examining the 
distributions of log(Ljv/L7v), where Lij is the central radio 
luminosity at the frequency of 1.4 GHz), Heckman (1995) 
argued that Sis are typically brighter than S2s by a factor 
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of ^4 (a result which is plagued by the large in- 
homogeneity of the radio fluxes); the distributions of 
log(LAr/L[o///]) yielded a lower factor of ~2, which is less 
inconsistent with ours. Our results are closer to those of 
Maiolino et al. (1995) and Giuricin et al. (1995), who 
derived a milder anisotropy from the direct comparisons 
of the Lpf luminosities. 

3.3 Relation between the IR properties and X-ray 
absorption 

As is expected within the framework of torus models, 

the NIR and MIR emissions of Scyfcrt nuclei appear to 
be related to their hard X-ray emission (see Danese et 
al. (1992), Kotilainen et al. (1992b) and Giuricin et al. 
(1995); Barcons et al. (1995), respectively). Furthermore, 
Goodrich et al. (1994) noted that the S2 galaxies in which 
broad Pa/3A1.2818/im emission lines have been detected 
are also the objects which generally have the lowest val- 
ues of the X-ray absorbing hydrogen columns Nh derived 
from the hard X-ray spectra. Recently, Granato et al. 
(1997) foimd some evidence that the X-ray column depths 
of S2 galaxies show a negative correlation with their N- 
band emissions, at least for objects which have a so low 



Compton depth {Nh <, 10 



24, 



that the hard X-rays 



can be seen directly through the torus. 

Our wide sample of IR data allows us to investigate 

further this point. To this end we have assembled the 
A^ij-values (mean values in the case of multiple entries) 
as derived by Smith & Done (1996) (for power-law fits 
to GINGA data), Iwasawa (1996), Warwick et al. (1993) 
for 17 Compton-thin S2 objects. Their A^/f-distribution 
ranges over 6 ■ 10^^ < Nh < 10^^ (with Nh in units of 
cm~^), with a median value of Nh = 7 • 10^^, which is 
considerably greater than the median value Nh ~ 10^^ of 
SI objects (Nandra & Pounds, 1994), even though the S2 
sample is biased towards low values of Nh- 

We examined the correlations between Nh and the 
IR colors discussed in §3.1 and those between Nh and 
the IR luminosities (normalized to the ^[0/77] luminosi- 
ties) discussed in §3.2. As regards the colors, only 
log(i^25/-Fi2) displays a significant (positive) correlation 
with Nh', the other colors show no correlations (but the 
number of objects is smaller). As regards the luminosities, 
log(Li2/£[0777]) shows a marginal anti-correlation with 
Nh (but log(L25/i[o777]) docs not for the same objects) 
and log(Li„/L[0777]) displays a strong anti-correlation; 
the lack of correlation for the other colors may simply be 
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due to the smaller number of objects. Fig. 4 illustrates the 
three above-mentioned significant correlations and Table 6 
presents the relevant results. In this table we list the two 
correlation coefficients, Spearman's r,, and Kendall's r^, 
together with the associated (one-tailed) percent signifi- 
cant levels (in these cases we have only detected objects). 

If we added in also the data for a few S2 galaxies which 
are likely to have Compton-thick {Nh > 10^'^cni~^) ob- 
scuring structures (i. e., NGC 1068, NGC 4945, NGC 
6552, and the Circinus galaxy; see, e.g.. Matt et al., 1996, 
and references cited therein) so that their hard X-ray 
emission is presumably due to a scattered nuclear com- 
ponent, all the correlations would vanish. This enlarged, 
less biased, sample of S2 objects has a median value of 
Nh ~ 2.9 • 1023cm-2. 

In summary, our study provides further support for the 
contention that the IR emitting torus is related to the X- 
ray absorbing structure, but only for Compton-thin struc- 
tures. 

4 Comparison with theoretical models 

The IR colors and the degree of anisotropy of the IR 

emission derived in the previous section can be compared 
with the relevant predictions of a wide series of torus mod- 
els (PK, GD, ER) characterized by different geometries, 
opacities and viewing angles. Wc have computed the the- 
oretical colors and anisotropics by filtering the predicted 
IR SEDs through the typical overall spectral responses 
(including filter transmission, detector response and at- 
mospheric transmission), which were taken from Bessell 
& Brett (1988) for the JHKL photometry, from Campins, 
Rieke & Lebofsky (1995) for the NQ photometry, from 
Beichman et al. (1988) for IRAS photometry. 

Models of compact cylinders 

PK modeled the torus simply as a geometrically thick, 
annular ring of uniform density dust which surrounds a 
central point source of radiation. PR proposed cylindri- 
cal configurations characterized by small extension (radial 
sizes <J a few pc) and very large optical thickness, in order 
to fit the general shape of the IR SED of Scyfcrt galax- 
ies (in particular that of the S2 galaxy NGC 1068) and 
the appearance of the ^ 10 /im silicate feature, which is 
seen in absorption in many S2 objects, whereas it is un- 
detected in many SI galaxies. Specifically, they preferred 
models with radial and vertical Thomson depths Tr=Tz=l 
(corresponding to a visual extinction Ay ^750 mag and a 
hydrogen column density A''// = l.S-lO^'^cm^^ for standard 
dust and dust-to-gas ratio) , an inner radius-to- height ratio 
a/h - 0.3, and T^s -500-1000 K (where T^s is the effec- 
tive temperature of the illuminating radiation on the inner 
wall of the torus). Basically, models with greater values 
of a/h or T^ff would shift the peak of the thermal bump 
to too short wavelengths and would lead to too warm IR 
colors, whereas models of low optical thickness (rr=0.1 
and/or Tz=0.1) would predict an (unobserved) prominent 
10 iJ,m emission feature for face-on viewing angles. 



Examining the predicted colors for face-on views (view- 
ing angle i = 0?), we have noticed that published SEDs (for 
Teff ~1000-2000 K) are in general characterized by too 
warm (small) IR colors, with SEDs which peak at shorter 
wavelengths than observed. Only models with low values 
of Te// {Teff ~500 K) and a/h {a/h ~0.1) can reasonably 
fit the MIR colors (i.e. log(F25/Fi2) and logiFg/ Fn)) of 
SI objects, for the wide range of values of Thomson op- 
tical depths (r ~0.1-10) explored by PK. But at shorter 
IR wavelengths those models imply a too steep IR SED, 
even for the choice of r-values which minimizes this prob- 
lem (i. e., Tz = Tr = 1); the relevant colors would require 
models with greater values of a/h and Te// {a/h —0.3-1, 
Teff ~1000 K). Similar problems hold for explaining both 
the long and short wavelength sides of the SEDs of S2s at 
the same time. 

Recognizing a problem which is emphasized by our wide 
data sample, PK noted that the \og{Fj\r / Flu) colors pre- 
dicted by their preferred models were appreciably greater 
than the observed colors of a small sample of 15 SI galax- 
ies. In order to make the observed mid-to-near IR slope of 
these galaxies and PC quasars fit better with a broadc;r IR 
bump, PR introduced a black body emission component 
ascribed to a small amount of hot (T ~1300 K) dust which 
lies inside the inner edge of the torus and is hidden by its 
main body from edge-on views. 

But the preferred optically thick tori, with Tz = Tr = 1, 
a/h = 0.3, Teff ~ 500 — 1000 K, entail an unreason- 
ably large degree of anisotropy in the IR emission. For 
instance, they predict a ratio between the luminosities re- 
ceived by a polar {i = 0°) and equatorial {i = 90?) observer 
L^{p)/L^{e) equal to ~4-7, ~10-20, -300-1500 at A ~25, 
lim , ~12 ^m , ~3.5 /xm , respectively. For cosi = 0.25 
in place of cost = 0, the corresponding ratios decrease re- 
spectively to ~2. 5-3.0, ~3-5, ~7-10, which are still large 
values. Only tori of low opacity, which however do not 
fit other characteristics, would predict a mild enough de- 
gree of anisotropy; for instance, tori with Tz — Tr = 0.1, 
a/h=0.3, Teff = 1000 K) imply L^{p)/L^{e) equal to ~1.7 
and ~2.5 at A ~ 25 /xm and ~12 /zm . 

4-2 Models of tapered disks 

In later models (GD, ER) inhomogeneous configurations 
with different outer shapes and density gradients were ex- 
plored. ER discussed three different types of extended 
models, i. e. flared discs, anisotropic spheres, and ta- 
pered discs. They soon realized that flared discs suf- 
fer the main drawback of predicting very prominent 10 
/xm emission features for face-on views. We have found 
that the two published models give too small colors for 
face-on views. For edge-on views, opaque models with 
equatorial optical depth Tuv = 200 (where Tuv = ta=ioooa) 
reproduce colors better than configurations with relatively 
low opacity (r„t,=40), except for log(F/v/Fi„), which still 
remains too low, but they predict too large an anisotropy 
of the IR emission (for standard dust and dust-to-gas ra- 
tio Av ^^0.24 Tuv)- The authors noted that very optically 
thick and compact flared disks would reduce the strength 
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of the emission feature, but they would also produce too 
narrow an IR continuum. 

Anisotropic spheres yield SEDs which show a more grad- 
ual transition from an edge-on to a face-on view than 
flared disks do. But face-on views still display strong 
10 /xm emission features and SEDs that are too flat in 
the MIR interval, whereas edge-on views, though pre- 
dicting fairly realistic colors for opaque configurations 
{tuv = 200), with the exception again of log(Fjv/-FLn) 
which remains too low, give too large an anisotropy. 

The authors devised tapered discs (discs whose height 
increases with distance from the central source but ta- 
pers off to a constant height in their outer part) mainly 
to suppress the 10 /xm emission silicate features from face- 
on views, without imposing an excessive equatorial depth 
(which would entail an unreasonably narrow IR contin- 
uum) or without resorting to lower silicate efficiencies (as 
proposed by Laor & Draine, 1993) or silicate depletion by 
shocks (as proposed by GD). The tapered disk models as- 
sume the following main parameters: the opening angle of 
the disk Oi, the ratio of inner to outer radii ri/r2, the ra- 
tio of the maximum half-thickness of the disc to the outer 
radius h/r2, the equatorial optical depth t„„, and the dust 
density distribution in radius. 

The authors asserted that configurations with large 
opening angles (8i ^ 60?) of the disk are unrealistic be- 
cause they tend to show 10 fim absorption features from 
all viewing angles, whereas configurations with small open- 
ing angles (9i ~ 30?) show implausible emission features 
for face-on views. Therefore, the authors deemed tapered 
disks with intermediate opening angles to be realistic and 
pubhshed the SEDs of four models with 6i = 45?, among 
which they preferred the one which predicts an almost fea- 
tureless MIR continuum together with a broad IR SED for 
face-on views. The preferred model is characterized by the 
following parametcrs:6i = 45?, t^v = 1500, density 
distribution, ri/r2=0.01, h/r2 = 0.3. 

We have found that none of the published models of 
tapered discs satisfactorily reproduce the observed small 
degree of anisotropy and the observed IR SEDs for face- 
on or edge-on views. In particular, the above-mentioned 
preferred model exhibits colors too large at longc;r IR wave- 
lengths and too small at shorter IR wavelengths for edge- 
on views, together with too large colors at shorter IR wave- 
lengths for face-on views. Besides, the model implies a 
large anisotropy corresponding to Lv{p) / Lv{e) equal to 
~10 and ^50 at A ~25 /xm and ~12 /zm . For a view- 
ing angle oi i = 60? the anisotropy, which remains large, 
decreases to ~2.5 and ~7 at A ~25 /xm and ~12 /xm , 
respectively. 

In conclusion, models of tapered disks, though having 
the virtue of solving problems related to the appearance 
of the 10 /tm feature, do not seem to easily reproduce the 
typical shape of the NIR and MIR SEDs of Seyfert nuclei 
and the observed anisotropy. 

Inspecting in detail NGC 1068, Efstathiou, Hough & 
Young (1996) realized that tapered disk models alone were 
unable to describe satisfactorily the IR SED of this object. 




1000 1500 

Fig. 5. — The plots show the colors predicted by torus mod- 
els as a function of the outer-to-inner radius rm/ro, for val- 
ues the equatorial optical thickness Te ranging from 20 to 80, 
in case of a face-on viewing angle. The two hatched regions 
indicate the 1 — <j and 2 — a error bars associated with the 
observed median. The left and right columns refer respec- 
tively to models with a standard silicate abundance and with 
a depression of silicates in the innermost regions. 
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Therefore, they supplemented the emission from the 
torus with an NIR component which was attributed to 
optically thin dust distributed in the ionization cone. In 
this way they obtained a more realistic SED, although the 
detailed fit of the 10 /um feature and of the extension of 
the NIR and MIR emissions are still difficult problems. 
We note that their composite model still implies a largely 
different SED for objects which are seen face-on and edge- 
on. 

4.3 Models of flared disks 

GD proposed flared disk conflgurations characterized by 
fairly large radial extension (from tens to hundreds pc), 
moderate optical thickness, an almost homogeneous den- 
sity distribution, i.e. tori with equatorial optical depth 
Te = TA3000A ~ 10 — 100 (for Standard dust and dust-to- 



gas ratio Ay ~ 0.61 • Te mag. and Nh ~ 1.2 • IO^^Te 



cm 



density distribution p oc with /3 < 0.5, outer-to- 
inner radius ratio Tm/ro ~ 300 — 1500 (where ro = 
0.5 (L/lO^^erg/s)^/'^ pc is set by the condition of subli- 
mation of all dust species). In this way GD were able to 
reproduce the average shape of the IR SEDs of a sample of 
16 optically-selected SI objects, without resorting to ad- 
ditional IR emission components. In order to suppress the 
occurrence of the 10 /^m silicate emission in these models 
for face-on views, GD proposed that shocks induced by 
radiation pressure depress the abundance of silicate grains 
in the innermost regions of tori (at r < few • lOro). 

Recently, from several arguments, which include the 
log FjM j Fj^n colors of a few nearby S2 galaxies and an 
upper limit on anisotropy in the MIR emission, Granato 
et al. (1997) argued in favor of moderately thick tori 
(with linc-of-sight extinction Ay <;80 mag) rather than 
of very thick configurations. Moreover, these authors fit- 
ted the IR SED of NGC 1068 with a line-of sight extinction 
Ay ^ 50 — 100 mag and with slightly less extended tori 
(with rm/ro ~ 100 — 200) than those proposed for Sis. If 
this were a general property of S2 galaxies, it would be a 
serious problem for the unified model of Seyfcrt galaxies. 
Clearly, this point deserves a statistical investigation such 
as the one we are carrying out in this paper. 

Tori of moderate optical thickness such as those pro- 
posed by GD certainly promise to be successful models for 
explaining our observational results, which lead to a quite 
low degree of IR emission anisotropy. Thus, we have de- 
cided to undertake a detailed comparison between a wide 
series of GD torus models and observational data, in or- 
der to assess to what extent these models arc compatible 
with our (partially new) results, which could provide more 
stringent constraints on torus properties. Using the GD 
code, we have computed the IR SEDs of many torus mod- 
els for which we have always assumed a uniform density 
distribution (/3=0) for a standard dust composition and a 
covering factor /=0.8, in order to keep the number of free 
model parameters to a minimum. We have parameterized 
this class of models by the two parameters rm/ro and Tg. 
The former parameter is related mainly to the broadness 
of the IR bump, whereas the latter controls the near-to- 



mid infrared slope of the SED along obscured lines of sight 
and the degree of its anisotropy. 

Fig. 5 shows the colors predicted by the GD torus mod- 
els as a function of rm/rg for different choices of the equa- 
torial optical thickness Te, in cases of face-on viewing an- 
gles. Fig. 6 show the analogous plots in cases of edge-on 
viewing angles. In Figs. 5 and 6 the two hatched regions 
mark the \-a and 2-a error bands associated with the 
observed median. We show models with a standard sili- 
cate abundance and models with a typical depression of 
silicate grains in the innermost regions (i.e., for r < 50ro). 
The latter models are intended to give an indication of the 
effects of a silicate depression on the colors. 




Fig. 6. — Same as in Fig. 5, but for an edge-on viewing 
angle. 

Let us first look at the models with a standard silicate 
abundance. They always tend to give values for the color 
log(F25/Fi2) that are greater than for \og{FQ / Fn) rather 
than vice-versa, as observed. The relatively high values of 
the latter quantity, namely the steepening of the SED as 
we go towards shorter wavelengths within the MIR spec- 
tral range, is not adequately reproduced by models. 
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Fig. 7. — For some photometric bands the plots show the 

ratio between the fluxes received by a polar and an equato- 
rial observer as a function of the equatorial optical thickness 
Te, for different values of the outer-to-inner radius rm/ro- 
The horizontal lines denote the median values derived from 
the comparison of the IR luminosities of SI and S2 nuclei. 
The left and right columns refer respectively to models with 
a standard silicate abundance and with a depression of sili- 
cates. 

It is thus reasonable to be content with reproducing 
a typical mean MIR color log(J25/-Fi2) ~ log(FQ/^jv) 
of ~0.3-0.4. These values rule out compact tori (i.e. 
with r„i/ro <^ 100). The color log(F/v/Fi„) is easily 
accounted for by models characterized by relatively low 
values of r^/rg and (e.g., r^/ro ~200 and Tg ~ 20) 
or by relatively high values of the two parameters (e.g., 
fm/fo ~700 and Tg ~80). In this case a decrease of rm/ro 
can be counterbalanced by an adequate decrease of Tg. A 
similar correlation between the two parameters holds for 
the color \og{FLn/ FKn), which rejects models with low 
rm/ro-vahica (e.g., with rm/ro ^300 for Tg ~20 or with 
rm/ro ^1000 for Te -80). The color \og{FKn/FHn) tends 
to reject tori of relatively small size and for very extended 
tori is not very sensitive to changes of rm/ro and Tg. The 
relatively low value of log(FH'„/Fj„) even more clearly fa- 
vors very extended tori {rm/ro ^ 500), within the usual 
Te-interval. 

The inclusion of silicate depression leaves the NIR SED 
almost unaffected, increases the MIR colors, and consid- 
erably reduces the color \og{F]y / F^n) ■ These effects are 
similar for all viewing angles. For face-on viewing angles 
the MIR colors are now in better agreement with the obser- 



vations. Now also log{F]^ / Fin) agrees with the quantities 
^og{FLn/FKn), ^og{FKn/ Fhu), log(i^//„/Fj„) in favoring 
very extended configurations. In summary, the IR colors 
better fit models of very extended tori(rm/ro ^500); Tg is 
not well constrained within the usual range (rg ~20-80). 
Notably, much larger and much lower values of Tg would 
give too steep and too flat MIR SEDs, respectively, for 
non-compact tori. 

Let us discuss edge-on systems, disregarding the unre- 
liable colors log(FKn/FHn) and log(i^H„/f j„). Tori of 
relatively small sizes {rm/ro ~ 100) can give realistic MIR 
colors and large enough values of \og{FM / FLn) for high 
values of Tg. Tori of intermediate sizes {rm/ro) do not fit 
the color log(i^7v/-pLri) very well. Very extended tori can 
again well fit all colors. The inclusion of silicate depres- 
sion does not change substantially the aforementioned con- 
siderations, although it tends to produce somewhat lower 
values of log(Fjv/.FL„) for very extended tori (for which, 
however, the assumption of uniform density distribution is 
likely to be less reasonable). But we have verified that the 
compatibility between the theoretical and observed values 
of log{FN / Fl,/) for very extended tori critically depends 
on the adopted degree of silicate depression. A lower de- 
gree of depression (i.e., a depression within r < 15ro) gives 
great enough values of that color. 

To Slim up, for S2s the two fit parameters are poorly 
constrained; the data clearly rule out only very compact 
tori {rm/ro ~10), which would require unreasonably large 
values of Tg. Within the unified schemes of AGNs, we 
would prefer a solution involving very extended tori also 
for S2 galaxies. Although the theoretical colors do not 
unambiguously converge towards this solution for S2s, it 
is encouraging that this possibility is not seriously contra- 
dicted. 

The comparison of the observed degree of anisotropy of 
the IR emission with the degree predicted by GD models 
provides a better constrain on Tg than that derived above. 
In Fig. 7 we show these predictions for some photomet- 
ric bands by plotting the ratio between the fluxes received 
by a polar and by an equatorial observer as a function of 
Tg, for different values of rm/ro- The horizontal lines de- 
note the median values (i.e. 1.5, 2.7, 6 for the 12 fim -, 
L-, K-bands, respectively) that we derived in §3.2 by com- 
paring the IR luminosities of SI and S2 nuclei. Models 
with a standard silicate abundance (left column) and with 
a standard depression of silicate grains (right column) in 
the innermost regions (for r < 50ro) yield similar results 
for moderate values of Tg. These plots favor values of Tg 
ranging from ~20 to —40 (i.e. 10 mag^ Ay ;^30 mag) for 
extended tori (with rm/ro ^ 100). 

5 Summary and Conclusions 

We have undertaken a statistical analysis of the MIR 
and NIR fluxes of an extended sample of Seyfert galaxies, 
the widest Seyfert sample so far considered in the relevant 
literature. We have paid particular attention in selecting 
homogeneous data (when available) and nuclear data, i.e 
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fluxes that arc, as far as possible, free of stellar light con- 
tamination. This has enabled us to derive the NIR and 
MIR continuum emission properties of Seyfert nuclei with 
unprecedented accuracy. We can summarize our main re- 
sults as follows. 

The color distributions indicate that the typical SED 
of S2 nuclei is always steeper than that of SI nuclei in 
the whole MIR and NIR spectral ranges. Hence S2s be- 
come typically fainter with respect to Sis as we go towards 
shorter IR wavelengths. 

There is no sure evidence of bimodality in the color and 
luminosity distributions. This agrees with the absence of 
significant bimodality in the distributions of several emis- 
sion properties of SI and S2 galaxies, such as the ultra- 
violet, soft X-ray, hard X-ray emissions (Mulchaey et al., 
1994) and radio core emission (Giuricin et al., 1996). A 
bimodality in some properties of S2s might be expected if 
there were two types of S2s, S2s with obscured BLR and 
without BLR, as was proposed by some authors (e. g., 
Hutchings & Ncff, 1991; Moran et al., 1992). 

The comparison of the IR luminosities of Sis and S2s, af- 
ter they arc made free of distance selection effects, reveals 
that S2s are somewhat brighter than Sis in the IRAS A ~ 
25 /um bands and in the IRAS FIR bands for the overall 
sample (though not for all subsamples). A direct compari- 
son of the IR luminosities, normalized to those of the [OIII] 
emission line, confirms these results. This finding, which 
is unexpected within the framework of the unified model 
of Seyfert galaxies, can be explained in terms of observa- 
tional selection effects favoring the inclusion of S2s with 
ongoing star formation in several current Seyfert samples, 
which are often dominated by optically-selected objects. 

As we go to shorter wavelengths, S2s become increas- 
ingly fainter than Sis by an amount which is in fairly 
good, quantitative agreement with what we expect from 
the differences between their SEDs. Within the frame- 
work of torus models, this behavior is consistent with an 
increasing anisotropy of a torus emission towards shorter 
IR wavelengths. But the degree of anisotropy we derive 
is quite low (i.e. lower than what has frequently been 
claimed in the literature). 

For a small sample of S2 galaxies having Compton-thin 
obscuring structures at hard X-ray energies, we discover 
significant correlations (previously unexplored) between 
the absorbing hydrogen columns derived from their hard 
X-ray spectra and some IR colors and luminosities. This 
indicates that the IR emitting torus is somehow related to 
the X-ray absorbing structure. 

The typical shapes of the IR SEDs of Sis and S2s and 
the luminosity differences between the two types of objects 
(i.e. the degree of anisotropy of the IR radiation) are com- 
pared with the predictions of a wide series of torus mod- 
els. Contrary to earlier claims, many models, in particular 
those involving compact or very thick configurations, arc 
found to be basically inconsistent with the observational 
data. 

We have specifically inspected the basic version of the 
torus models by Granato & Danese (1994), which was in- 



tended to represent the simplest kind of tori envisaged by 
the AGN unification schemes. We have found that the ma- 
jor IR observational properties can be reasonably satisfied 
by flared disk configurations characterized by quite a large 
extension (a few hundreds pc) and moderate equatorial op- 
tical thickness (corresponding to a typical extinction of 10 
mag^ Ay ^30 mag), for both SI and S2 galaxies. 

For a standard dust over gas ratio, the aforementioned 
range of torus optical thickness, derived from fitting the 
IR properties of Seyfert nuclei, translates into lower val- 
ues of Nh (Nh ~ 2 • 10^^ - 5 • lO^^cm-^) than the median 
value of A''^ ~ 2.9 • lO^^cm"^ of a sample of S2 galax- 
ies, which is still probably biased towards low values of 
Nh- This lends further support to the recent view that 
the matter responsible for the X-ray absorption is mostly 
dust-free and that the bulk of the X-ray absorption does 
not occur in the dusty torus, but at and just inside the 
dust sublimation radius (Granato et al., 1997). 

Extended tori arc consistent with the detection of dusty 
nuclear disks on scales of '^lOO pc in some Hubble Space 
Telescope observations (e.g., Jaffe et al., 1993) and with 
the detection of dense molecular nuclear disks on simi- 
lar scales in some high resolution radio observations (e. 
g., Tacconi et al., 1994; Kohno et al., 1996). Extended 
tori are also consistent with the presence of red (i.e. with 
V — R > 1.0 mag) regions of hundreds of pc which co- 
incide with the optical nuclei of many S2 objects; they 
can be interpreted as being due to reddening by dusty tori 
(Mulchaey et al., 1996b). It is known that there is a se- 
lection bias against finding AGNs (especially UV-excess 
AGNs and spectroscopically-selected AGNs) in highly in- 
clined galaxies. This has led some authors to suggest that 
the obscuring material covers a part of the NLR (over 
;^100 pc) and is roughly coplanar with the host galaxy 
disk (McLeod & Rieke, 1995; Maiohno & Rieke, 1995). 

Extended configurations are also in better agreement 
than compact tori (with a radial size of ;^10pc) with the 
high-resolution NIR and MIR observations of NGC 1068, 
which showed an emission extended over '-^l" scale (e.g., 
McCarthy et al., 1982; CheUi et al., 1987; Tresch-Fienberg 
et al., 1987; Braatz et al., 1993; Cameron et al., 1993). If 
this emission were due to torus radiation reflected by elec- 
tron scattering and to NLR dust (Pier & Krolik, 1993), 
there would probably be too large reddening in the NLR 
(Cameron et al., 1993). The analysis of NIR polarized 
images of the nuclear region of this galaxy points to a 
torus with a size of ~200 pc (Young et al., 1996b). An 
extended (~2") bar- like NIR emission, which might hint 
at an extended dusty torus viewed approximately edge- 
on, has been recently detected in Mk 348 (Simpson et al., 

1996) . 

Tori with moderate optical depths are consistent with 
the detection of broad components of NIR lines in sev- 
eral narrow-line AGNs (as discussed by Granato et al., 

1997) . In most cases of detections of broad lines in a seem- 
ing narrow- line object, the broad lines are reddened more 
than the narrow lines, locating the bulk of the dust re- 
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sponsible for absorbing the broad-linc emission between 
the BLR and the NLR (e.g., Hill, Goodrich & DePoy, 
1996). Narrow-line galaxies which do not show broad 
lines in total or polarized flux could have obscured regions 
of scattered radiation or lines which are thermally very 
broadened and hence hardly observable (see, e.g.. Young 
et al., 1996a). Moderate optical depths are also consistent 
with the typical extinctions (Ay 30 mag) inferred from 
the strengths of the 10 /xm silicate absorption features ob- 
served in several S2 galaxies (e.g., Aitken & Roche, 1985; 
Roche et al., 1991). 

Alleviating some earlier worries, we can conclude that 
the general unified scheme of AGNs is basically compatible 
with the IR observational data. Notably, this scheme can 
not hold in a strict sense, since the broadness of the color 
distributions suggests a large spread in the properties of 
obscuring tori, especially for S2s. Other arguments against 
the strict version of the unified scheme can be found in An- 
tomicci's (1993) review. 

There is no indication in our sample that the optical 
depth of tori decreases with AGN luminosity, as is ad- 
vocated by the simplest modification of the strict unified 
scheme (e.g., Lawrence, 1991). By the way, the SEDs of 
some hyperluminous high-redshift IR galaxies have been 
recently reproduced in terms of a luminous AGN sur- 
rounded by an extended torus with a larger equatorial 
optical depth (60 mag ^ Ay <^ 300 mag) than those of 
classical Scyfcrt galaxies (Granato, Danese & Franceschini, 
1996). The ISO satellite observations are expected to sub- 



stantially improve our knowledge of the IR properties of 
high-luminosity AGNs. 

Nevertheless, our detailed comparison with the obser- 
vational data raises several problems. Also the most suc- 
cessful torus models, i.e. extended and moderately thick 
flared disks, fail to reproduce accurately several observa- 
tional details. Even with the present quality and quantity 
of observations there is a need for more sophisticated ver- 
sions of torus models having more free parameters to fit the 
data (e.g., different matter distributions and geometries), 
though not so many as to preclude physical understand- 
ing. Remarkably, especially in the case of extended con- 
figurations, the possible clumpy nature of the actual dust 
distribution weakens the validity of current torus models, 
which are based on a continuous distribution of matter. 



The authors are grateful to L. Danese for enlightening 
discussions on torus models and to M. Mezzetti for com- 
ments on the statistical analysis. 

The authors thank E. D. Feigelson for having kindly 
provided the ASURV software package. 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED), which is operated by the Jet 
Propulsion Laboratory, Caltech, under contract with the 
National Aeronautics and Space Administration. 

This work has been partially supported by the Italian 
Ministry of University, Scientific, and Technological Re- 
search (MURST) and by the Italian Space Agency (ASI). 



REFERENCES 



Aitken, D. K. & Roche, P. F., 1985, MNRAS, 213, 777. 
Alonso-Herrero, A., Ward, M. J. & Kotilainen, J. K., 1996, MNRAS, 

278, 902. 

Antonucci, R. R. ,T., ARA&A, 31, 473. 

Barcons, X., Franceschini, A., De Zotti, G., Danese, L. & Miyaji, T., 

1995, ApJ, 455, 480. 
Barvainis, R., 1987, ApJ, 320, 537. 
Barvainis, R., 1990, ApJ, 353, 419. 

Barvainis, R., 1992, in "Testing the AGN paradigm", S. Holt, S. Neflf 

& M. Urry eds. (New York: AIP), 129. 
Beichman, C. A., Neugcbaucr, G., Habing, H. J., Clegg, P. E. & 

Chester, T. J., 1988, IRAS Catalogs and Atlases: Explanatory 

Supplement (Washington, DC: GPO). 
Bessell, M. S. & Brett, J. M., 1988, PASP, 100, 1134. 
Binette, L., Fosbury, R. A. & Parker, D., 1993, PASP, 105, 1150. 
Blanco, P. R., Ward, M. J. & Wright, C. S., 1990, MNRAS, 242, 4P. 
Boisson, C. & Durret, F., 1986, A&A, 168, 32. 

Bothun, G. D., Halpern, J. P., Lonsdale, C. J., Impey, C., Schmitz, 

M., 1989, ApJS, 70, 271. 
Braatz, J. A., Wilson, A. S., Gezari, D. Y., Varosi, F. & Beichman, 

C. A., ApJ, 409, L5. 
Cameron, M. et al., 1993, ApJ, 419, 136. 

Campins, H., Rieke, G. H. & Lebofsky, M. J., 1985, AJ, 90, 896. 
Carico, D. P. et al., 1988, AJ, 95, 356. 

ChcUi, A., Perricr, C, Cruz-Gonzales, I. & Carrasco, L., 1987, A&A, 
177, 52. 

Dahari, O. & Dc Robertis, 1988, ApJS, 67, 249. 

Danese, L., Zitelli, V., Granato, G., Wade, R., De Zotti, G. & Man- 

dolesi, N., 1992, ApJ, 399, 38. 
de Grijp, M. H. K., Keel, W. C, Miley, G. K., Goudfrooj, P. & Lub, 

J., A&AS, 96, 389. 
Edelson, R. & Malkan, M. A., 1987, ApJ, 323, 516. 
Edelson, R., Malkan, M. A. & Rieke, G. H., 1987, ApJ, 321, 233. 



Efstathiou, A., Hough, J. H. & Young, S., 1996, MNRAS, 277, 1134. 
Efstathiou, A. & Rowan-Robinson, M., 1995, MNRAS, 273, 649 
(ER). 

Elvis, M., Willner, S. P., Fabbiano, G. et al., 1984, ApJ, 280, 574. 

Evans, I. N., Koratkaj", A. P., Storchi-Bergmann, T. et al., 1996, 
ApJS, 105, 93. 

Feigelson, E. D. & Nelson, P. I., 1985, ApJ, 293, 192. 

Gezari, D. Y., Schmitz, M., Pitts, P. S. & Mead, J. M., 1993, "Cat- 
alog of Infrared Observations" , NASA Reference Publication No. 
1294. 

Giuricin, G., Fadda, D. & Mezzetti, ApJ, 468, 475. 

Giuricin, G., Mardirossian, F. & Mezzetti, M., 1995, ApJ, 446, 550. 

Giuricin, G., Tamburini, L., Mardirossian, F., Mezzetti, M. & 

Monaco, P., 1994, ApJ, 427, 202. 
Glass, I. S., 1992, MNRAS, 256, 23 p. 

Glass, I. S. & Moorwoood, A. F. M., 1985, MNRAS, 214, 429. 
Golombek, D., Miley, G. K. & Neugebauer, G., 1988, A J, 95, 26. 
Goodrich, R. W., Veilleux, S. & Hill, G. J., 1994, ApJ, 422, 521. 
Granato, G. & Danese, L., 1994, MNRAS, 268, 235 (GD). 
Granato, G., Danese, L. & Franceschini, A., 1996, ApJ, 460, Lll. 
Granato, G., Danese, L. & Franceschini, A., 1997, ApJ, in press. 
Heckman, T. M., 1995, ApJ, 446, 101. 

Heisler, C. A., De Robertis, M. M. & Nadeau, D., 1996, MNRAS, 
280, 579. 

Helou, G., Khan, I. R., Malck, L. & Boehmcr, L., 1988, ApJS, 68, 
151. 

Hes, R., Barthel, P. D. & Fosbury, R. A., 1993, Nature, 362, 326. 
Hill, G. J., BeckUn, E. E. & Wynn-Williams, C. G., 1988, ApJ, 330, 
737. 

Hill, G. J., Goodrich, R. W. & DePoy, D. L., 1996, ApJ, 462, 163. 
Hoel, P. G., 1971, Introduction to Mathematical Statistics (4th ed.: 

New York, Wiley). 
Hutchings, J. B. & Neff, S. G., 1991, AJ, 101, 434. 



17 



Isobe, T., Feigelson, E. D. & Nelson, P. J., 1986, ApJ, 306, 490. 

Iwasawa, K., 1996, Ph.D. thesis, ISAS. 

Jackson, N. & Browne, I. W., 1990, Nature, 343, 43. 

Jaflfe, W., Ford, H. C, Fcrraresc, L., Van den Bosch, F. & O'Connell, 

R. W., 1993, Nature, 364, 213. 
Laing, R. A., Jenkins, C. R., Wall, J. V. & Ungcr, S. W., 1994, in 

ASP Conf. Scr. 54, The First Stromlo Symposium: The Physics 

of Active Galaxies, eds. G.V. Bicknell, M. A. Dopita & P. J. 

Quinn (San Francisco: ASP), 175. 
Lawrence, A., 1991, MNRAS, 252, 586. 

Keel, W. C, de Grijp, M. H. K., Miley, G. K. & Zheng, W., 1994, 
A&A, 283, 791. 

Kendall, M. & Stuart, A., 1977, The Advanced Theory of Statistics 

(London: Griffin). 
Knapp, G. R., Bies, W. E. & van Gorkom, J. H., 1990, AJ, 99, 476. 
Knapp, G. R., Guhathakurta, P., Kim, D.-W. & Jura, M., 1989, 

ApJS, 70, 329. 

Kohno, K., Kawabe, R., Tosaki, T. & Okumura, S. K., 1996, ApJ, 
461, L29. 

Kotilainen, J. K. & Prieto, M. A., 1995, A&A, 295, 646. 
Kotilainen, J. K., Ward, M. J., Boisson, C., DcPoy, D. L., Bryant, 

L. R., Smith, M. G., 1992a, MNRAS, 256, 125. 
Kotilainen, J. K., Ward, M. J., Boisson, C., DePoy, D. L., Bryant, 

L. R., Smith, M. G., 1992b, MNRAS, 256, 149. 
Laor, A. & Drainc, B. T., 1993, ApJ, 402, 441. 
La Valley, M. Isobe, T. & Feigelson, E. D., 1992, BAAS, 24, 839. 
Lawrence, A. et al., 1985, ApJ, 291, 117. 

Loska, Z., Szczerba, R. & Czcrny, B., 1993, MNRAS, 261, 63. 
Maiolino, R. & Rieke, G. H., 1995, ApJ, 454, 95. 
Maiolino, R., Ruiz, M., Rieke, G. H. & Keller, L. D., 1995, ApJ, 446, 
561. 

Matt, G., Fiore, F., Perola, G. C. et al., 1996, MNRAS, 281, L69. 
Mazzarella, J. M., Bothun, G. D. & Boroson, T. A., 1991, AJ, 101, 
2034. 

McAlary, C. W., McLaren, R. A. & Crabtree, D. R., 1979, ApJ, 234 
471. 

McAlary, C. W., McLaren, R. A., McGonegal, R. J., Maza, J., 1983, 

ApJS, 52, 341. 
McAlary, C. W. & Rieke, G. H., 1988, ApJ, 333, 1. 
McCarthy, D. W., Low, F. J., Kleinmann, S. G. & Gillet, F. C., 

1982, ApJ, 257, L7. 
McLeod, K. K. & Rieke, G. H., 1995, ApJ, 441, 96. 
Miles, J. W., Houck, J. R., & Hayward, T. L., 1994, ApJ, 425, L37 
Moran, E. C., Halpcrn, J. P., Bothun, G. D. & Becker, R. H., 1992, 

AJ, 104, 990. 

Moshir, M. et al., 1992, Explanatory Supplement to the IRAS Faint 

Source Survey, Ver. 2 (Pasadena: JPL). 
Mulchacy, J. S., Koratkar, A., Ward, J. et al., 1994, ApJ, 436, 586. 
Mulchacy, J. S., Wilson, A. S. & Tsvetanov, Z., 1996a, ApJS, 102, 

309. 

Mulchaey, J. S., Wilson, A. S. & Tsvetanov, Z., 1996b, ApJ, 467, 
197. 

Nandra, K. & Pounds, K. A., 1994, MNRAS, 268, 405. 
Nelson, C. H. & Whittle, M., 1995, ApJS, 99, 67. 
Neugebauer, G., Green, R. F., Matthews, K. et al., 1987, ApJS, 63, 
615. 

Neugebauer, G., Soifer, B. T., Matthews, K. & Ehas, J. H., 1989, 
AJ, 97, 957. 

Osterbrock, D. E. & Martel, A., 1993, ApJ, 414, 552. 

Phinney, E. S., 1989, in Theory of Accretion Disks, ed. J. S. Meyer 



(Dordrecht: Kluwer), 457. 
Pier, E. A. & Krolik, J. H., 1992, ApJ, 401, 99 (PK). 
Pier, E. A. k, Krolik, J. H., 1993, ApJ, 418, 673. 
Pogge, R. W., 1989, ApJ, 345, 730. 

Rice, W., Lonsdale, C. J., Soifer, B. T. et al., 1988, ApJS, 68, 151. 
Rickc, G. H., 1978, ApJ, 226, 550. 

Roche, P. F., Aitken, D. A., Smith, C. H. & Ward, M., 1991, MN- 
RAS, 248, 606. 

Rudy, R. J., Levan, P. D. & Rodriguez-Espinosa, J. M., 1982, A J, 
87, 599. 

Rudy, R. J. & Rodriguez-Espinosa, J. M., 1985, ApJ, 298, 614. 
Ruiz, M., Rieke, G. H. & Schmidt, G. D., 1994, ApJ, 423, 608. 
Rush, B., Malkan, M. A. & Spinoglio, L., 1993, ApJS, 89, 1. 
Salzer, J. J. & MacAlpine, G. M., 1988, AJ, 96, 1192. 
Sanders, D.B., Egami, E., Mirabel, I. F. & Soifer, B. T., 1995, AJ, 
110, 1993. 

Sanders, D. B., Phinney, E. S., Neugebauer, G., Soifer, B. T. & 

Matthews, K., 1989, ApJ, 347, 29. 
Shapiro, S. S. & Wilk, M. B., 1965, Biometrika, 52, 591. 
Simpson, C., Mulchaey, J. S., Wilson, A. S., Ward, M. J. & Alonso- 

Herrero, A., 1996, ApJ, 457, L19. 
Smith, D. A. & Done, C., 1996, MNRAS, 280, 355. 
Soifer, B. T., Boehmer, L., Neugebauer, G. & Sanders, D. B., 1989, 

AJ, 98, 766. 

Spinoglio, L. & Malkan, M. A., 1989, ApJ, 342, 83. 

Spinoglio, L., Malkan, M. A., Rush, B., Carrasco, L. & Recillas— 

Cruz, E., 1995, ApJ, 453, 616. 
Stein, W. A. & Wcedman, D.W., 1976, ApJ, 205, 44. 
Tacconi, L. J., Genzel, R., Blictz, M. et al., 1994, ApJ, 426, L77. 
Telesco, C. M., 1988, ARA&A, 26, 343 
Telesco, C. M., & Decher, R., 1988, ApJ, 334, 573 
Telesco, C. M., Dressel, L. L., & Wolstencroft, R. D., 1993, ApJ, 

414, 120 

Trcsch-Ficnbcrg, R. et al., 1987, ApJ, 312, 542. 
TuUy, R. B., 1988, Nearby Galaxies Catalog (Cambridge: Cambridge 
Univ. Press). 

Vader, P. J., Progel, J. A., Terndrup, D. M. & Heisler, C. A., 1993, 
AJ, 106, 1743. 

Veron-Cetty, M. P. & Veron, P., 1991, A Catalogue of Quasars and 

Active Nuclei, 7th ed. (ESQ Sci. Rep. No. 17). 
Ward, M. et al., 1987, ApJ, 315, 74. 

Warwick, R. S., Sembay, S., Yaqoob, T. et al., 1993, MNRAS, 265, 
412. 

Whittle, M., 1992, ApJS, 79, 49. 
Willner, S. P. et al., 1984, PASP, 96, 143. 

Wilson, A. S., 1992, in "Physics of Active Galactic Nuclei", eds. W. 
J. Duschl & S. J. Wagner (Springer Verlag: Heidelberg), 307. 

Wilson, A. S., 1996, Vistas in Astronomy, 40, 63. 

Wilson, W. J., Schwartz, P. R., Neugebauer, G., Harvey, P. M. & 
Bccklin, E. E., 1972, ApJ, 177, 523. 

Wynn- Williams, C. G. & Becklin, E. E., 1993, ApJ, 412, 535. 

Young, S., Hough, J. H., Efstathiou, A., Wills, B. J., Bailey, J. A., 
Ward, M. J. & Axon, D. J., 1996a, MNRAS, 281, 1206. 

Young, S., Packham, C, Hough, J. H. & Efstathiou, A., 1996b, MN- 
RAS, 283, LI. 

Zhou, S., Wynn-WiUiams, C. G. & Sanders, D. B., 1993, ApJ, 409, 
149. 

ZiteUi, v., Granato, G., Mandolesi, N., Wade, R. & Danese, L., 1993, 
ApJS, 84, 185. 



18 



